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Charge kinks are considered as fundamental excitations in quarter-filled charge-ordered ladders. 

The strength of the coupling of the kinks to the three-dimensional lattice depends on their energy. 

The integrated intensity of Raman scattering by kink-antikink pairs is proportional to (f>^ or 0'*, 

where is the order parameter. The exponent is determined by the system parameters and by the 

strength of the electron-phonon coupling. 
An internlay of spin, charge, and lattice degrees of freedom drives a variety of phase transitions in ladder-like 
compoundstl. The best known examples of these systems are half-filled SrCu203 and Sri4Cu2404i, where the mag- 
netic properties are determined by s = 1/2 spins at the Cu sites with a superexchange interaction arising from 
hopping via oxygen orbitals. Ca-substitution in Sr sites introduces mobile holes and leads to superconductivity in 
Sr2.5Caii.5Cu2404i under pressureH. The pairing interaction has a magnetic origin since it isi favorable to build up a 
local spin singlet of two holes if the exchange interaction overcomes the kinetic energy effectgj. In Sri4_a;Caa;Cu2404i 
the dopant holes reside at the oxygen rather than at the Cu orbitals responsible for magnetic properties. Another 
class of ladders are the half-filled vanadates such aSpCaV205 and MgV205, where spins 1/2 are located at the V sites. 
A more interesting realization is given by NaV205a and LiV205, which have VO5 pyramids as a common structure 
element with the half-filled vanadates. As follows from the chemical formula, the mean charge of V in NaV205 and 
LiV205 is 4.5. Although similar in stoichiometry and unit cell structure, these two compounds have strikingly different 
physical properties. In LiV2 05, the V ions are arranged in parallel chains of non-magnetic and magnetic V^^. In 
NaV205 all V ioiiS|-are equivalent, making it a rare realization of the quarter-filled laddeiolj. Only for temperatures 
below Tco ~ 34 KB'lI the compound undergoes a charge ordering transition into a spin-gapped phase, which is due to 
an interplay of charge, spin, and lattice degrees of freedftaaE. It has been proposed that the spin-lattice interaction in 
NaV2 05EI is due to a strong electron-phonon couplingOU. 

Low-energy excitations in the ladders related to charge and spin degtees. of freedom and having energies less than 
0.2 eV, are well suited for investigations by light scattering spectroscopyHEj. However, the mechanism of their Raman 
activity, being strongly related to the way they modulate the crystal's polarizability, requires a detailed investigation 
for each excitation. 

Since the ladders are one-dimensional systems with strong-i pier actions, non- linear models appeared to be very 
useful for investigations of charge- and lattice-llJ, and magneticollj excitations. At small doping, when the interaction 
between the carriers is npt very strong, the generalized t — J model already contains the essential physics necessary 
to describe the ladderscl. In view of the large doping in the quarter filled ladders (a depletion of 1/4 hole per 
magnetically active site) the excitations in the charge ordered state should again be considered within a non-linear 
framework. In this paper, we study the excitations in a charge ordered quarter-filled ladder system with strong 
coupling to the lattice. The elementary excitations are charge kinks involving ion displacements. We establish a 
mechanism for Raman scattering from these excitations and calculate the scattering intensity which is found to be 
strongly dependent on the order parameter. 

In a quarter-filled ladder each rung i is occupied by a charge Qi — (5,|^^ + Q-^^ = 1, where \L) and \R) denotes the 
left and right leg state, respectively. The electronic degree of freedom is the hopping of the charge between two ions 
of a rung described by a matrix element t±_ . Hopping between nearest-neighbor rungs is prohibited since it requires 
the energy of approximately 2t±^ to put electrons from two binding two-site states, with the energies of —t^ each, to 
two one-site states at one rung. The wave function of an electron on a rung has the form ipi — aL\L) -\- aii\R), where 
a\ + a\ = 1. In the disordered phase = —cur = l/\/2- 

The total energy of the system Etot = t^cor + K + U\at is the sum of the correlation energy J7cor due to the 
repulsion of electrons at nearest-neighbor sites described by V^or, the kinetic energy and the lattice term C/iat 
which includes the electron-phonon coupling. Ucot is the sum over the nearest neighbor interactions of the type 

Kor ( Q^^QI+i + ) ■ With 0i = 0f = Qf - 1/2, this contribution can be written as: /7cor = Kor I]^ ^i^i+i, 
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where (f>f = -=jif has been used. The kinetic energy per rung can be expressed as: Ki = —t± + 2t± (<?!)? + (j)^^ . The 
lattice terrnQ't2l has the Holstein-like form Uut = C'^cj)[zij + The elastic contribution to the energy due 

to a displacement Zij of the ions at rung i and leg I is 1/2 ; '^^ih with lattice force constant k = Mft^, where M 
is the mass of the atom, f2 the vibrational frequency, and C is the deformation potential. The static ion displacement 
minimizing [/lat is Zi = - (C / k) (jji. n ITl 

The equilibrium state of the system minimizes -Etot- The result is a zigzag ordered statetZI"E2l with the order 
parameter = 0iat ■ (—1)% where cj)^^^ = (2(T4or — ^.l) + Vut) /4:t±, and Mat = C^/k. The subscript "lat " serves as 
a reminder that the lattice is involved in the charge ordering. The phase transition occurs when 2Vi3or + Vy^t > '2t±. 

In order to study the dynamics, the energy i^tot is rewritten in a continuum approximation, with a time-dependent 
order pMameter 4>{y,t) that slowly changes along the ladders. The energy acquires a form typical for the standard 
04 fieldEl: 



E\ 



{d^/sutdtf + {dcj^/dyf + A _ ^l^f /2] ^, (1) 



where l/sf^t = MC^/KorolK^ with ay being the lattice constant along the ladder, and A — 4i:j_/Koray • The 
d(j)/dt term in Eq.(l) arises due to the kinetic energy of the V ions which displacements follow their charges since 
dzi/dt = — {G / k) d(j>i/dt. The ground state for E [(f)] is 4>{y) = ±0iat. The classical excitations within this model are 
kinks and antikinks of the form: 

Hy^t) = ±(/)iattanh0iat\/ • —^==^=. (2) 

V 2 ^1 _ uysi. 

Here u is the kink velocity, and its energy Eu = 7iat-Eiat, where E^t = 4:y/2^/t± V~(j)f^j3, and 7iat = 1/^1 - w^/sfaf 
The density of states of the one-kink excitation is given by: p{E) = E/T:sut\/ E^ — Ef^^. The displacements of ions 
in the ground state and in a kink excitation are shown in Fig. 1. 

The above consideration is valid for a "soft" lattice, where ion displacements follow the electron redistribution 
during the kink propagation. The "soft" lattice condition is r^ri 3> 1, where the time t\ characterizes the rate of 
the change of the order parameter. It is determined by the kink width Wk ^ l/Tiat^iatV^ as Tk ^ w-k/u. In the 
"ultrarclativistic" limit, where u siat, the condition Tk^l ^ 1 cannot be fulfilled and the lattice becomes "rigid" 
due to the Lorentz contraction of the kink width ^ 1 /7iat . 

To understand the difference between a soft and a rigid lattice, let us consider the condition Tkil 3> 1 in more detail. 
Let us assume first that the "light" velocity siat is large enough so that fJ/^iatA/Asiat «C 1. In this case the limiting 
velocity Umax, determined by the condition Tk ^ i^"^, is Umax ~ ^^/0iat%/A <C Siat- At u > Umax the kinks become 
decoupled from the lattice and propagate on the background of ions displaced as in the equilibrium charge ordered 
state. The order parameter for the rigid lattice — (V^or ~ t±)/'^t± is determined by the electronic subsystem only, 
and the "light" velocity Sd 3> siat since the ions are not involved in the kink motion anymore. It seems that the 
excitation energy drops down to approximately E^i = £'iat0ei/</'fat ^^"^ becomes weakly u— dependent. However, due 
to the increase of the "light" velocity, the Lorentz contraction disappears and the soft lattice regime is restored. This 
fact implies that the lattice becomes rigid only for the fast kinks with u close to fl/(f>eiV^- In other words, there are 
no kinks in the interval il/0iat%/A J; u < ri/0ci\/A. 

In the opposite case when r2/0iat\/Asiat 2> 1, the lattice becomes rigid when Umax is very close to siat, such that 

-^lat Umax ^ -^lat ^ri/0iatVAsiat^ • Therefore, the kinks involving lattice displacements are well defined up to high 

energies Eu^^^ ^ £'iati^/</>iatVAsiat. Here the gap in the allowed u is determined by siat < u < Vt/(t>c\\/\. If K;or < 
the kink-like electronic excitations are not well defined being strongly overdamped, and the "rigid lattice" electronic 
part of the spectrum is absent. 

Now consider light scattering by the kinks. These excitations modulate the charge density and the crystal's dielectric 
function e^^'', where w is the light frequency, and (i and rj are Cartesian indices, thereby causing inelastic light scattering 
at frequencies equal to the excitation energy. 

The variation of the dielectric function is proportional to the square of the order parameter, and can be written for 
one rung as (the Cartesian indices are omitted): 

m ~ (0L) = ^ (0? - '^fat) • (3) 
The change of the polarizability per kink 4>{y,t) from Eq.(2) is: 
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Ae.„. = ^f my,t) - ^ = -2V2||^. (4) 

J n "'^ TlatVA 

With increase of the energy, Aekink decreases due to Lorentz contraction of the kink width. 

Because of the very small photon wavevector, Raman scattering probes excitations with zero net momentum. 
Therefore, the Raman active quasiparticles are kink-antikink pairs with velocities u and —u, respectively. The 
corresponding contribution of the kink-antikink pairs to the polarizability e^** is shown in Fig. 2 for two different kink 
energies. The measured spectral density of the scattered light p{E) as a function of energy transfer to the system 
(Raman shift) is proportional to the probability of the excitation of the pair with the total energy E = 2i?„. Since 
the Raman scattering is a process of the decay of the incident photon into a continuum consisting of the scattered 
photons and electronic excitations, its probability is given by Fermi's Golden Rule as: 

Pkink(2i?) = 27rK£;) (Aekinkf = 16 f %^=- (5) 

\drj Asiat7iat\/7iit-l 

The spectral density of the scattered light in Eq.(5) has a threshold at 2i?iat0 and decreases at i? 3> E\^t as E^^ . 
The integrated intensity of the Raman continuum of kinks in the "soft" lattice can be written as: 



T71 • ' ""^ 



/kink = / pi2E)dE = Stt — ^ v^^iat. (6) 



ASl; 



The integration in Eq.(6) was extended to infinity leading to /kink '/'fat/'iat ~ ^fat- This can be done if 
/^lat "v/Asiat ^ 1, which is, as we will see below, the case for NaV205. In the opposite case of ri/0iatA/Asiat 1, the 
integration should be performed up to E ^f2/0iat\/A^ , which is close to /?iat that yields /kink ~ '/'fat- 

Let us now-discuss a contribution of the kink-antikink excitations to the Raman spectra of NaV205, which contains 
a continuunJiju and several phonon peaks. The changes at T < Tco manifest themselves by changes in the phonon 
frequencies and line shapes, and by the appearance o£_tiew intense peaks. Some of these peaks are vibrational modes 
while others can be attributed to magnetic excitationalj. The modification of the continuum consists of two effects: a 
depletion at the spectral range up to 30 meV and a moderate redistribution of the spectral weight at higher energies. 
This observation implies that the Raman scattering mechanism in the ordered and disordered phases must be virtually 
the same, despite the different character of the excitations. 

The relevant part of the unit cell of the quarter filled ladder compound NaV2 05 is shown in Fig.3|. The Holstein-like 
electron-phonon coupling is due to the 03 ions located either above or below the V-O-V rungsQ with the hopping 
matrix element t± « 0.35 eV. Coupling to the lattice favors the charge ordering since l^at > 0, hcwever it does not 
drive the transition itself since in the absence of correlations the lattice is stable, that is Mat < 2t^a. 

The charge ordering parameter 0iat determined froHi|the magnetic susceptibility is approximately 0.20. Another 
possibility to find 0iat is provided by structural dataE3 where one finds l^il « 0.04 A. Assuming C = 5 eV/A, and 
n = 400 cm~^ one obtains k ~ 40 eV/A^, V|at « 0.6 eV, and (?!)iat = (n/C) Zi « 0.3. We accept this magnitude of 
(/)iat for further estimates and for t±_ = 0.35 eV, and Vcov = t± obtain E\at ~ 20 meV. Therefore, one can expect a 
depletion in the scattering intensity for the Raman shift less than 2E\at = 40 meV, in agreement with the experimental 
datall3. The chosen model parameters and a|| = 4A yield ri/(/iiat\/Asiat ~ 10, that is the lattice is soft up to high kink 
energies. This fact justifies an extension of the integration in Eq.(6) to infinity. Considering the kinks, we neglected 
the interaction between different ladders. As we established above, at Vcor < t±, the "rigid lattice" part of the kink 
spectrum is absent. We assume that this is the case in NaV205 since Vcor + Mat/2 — i± should be considerably smaller 
than tj_ to ensure a small order parameter. For this reason we did not consider scattering by the high-energy kinks 
uncoupled to the lattice. 

To have a reference point, we compare the spectral densities of the two-kink scattering and two other Raman 
processes relevant for this compound. The first one is the first-order phonon Raman scattering. The estimate of 
spectral density at the phonon peak maximum is given by: pph ~ (de^/dQ) {zq/cl) a||/r, where Q is the ratio of the 
ion displacement to the characteristic lattice constant a, zq is the zero-point vibrational amplitude for the mode, and F 
is the phonon linewidth. The other mechanism is two-magnon Raman scattering due to the frustrated superexchange 
interaction. Although the scattering intensity depends on many detailfj, an estimation within the Fleury-Loudon 
theoryn3 can be done as: Pmagn ^ {J/Eo)'^a^^/J, where J is the superexchange and Eq is of the order of magnitude of 
the interband transition energy. By taking into account that siatOn"^ ^ Q, the ratio of intensities can be estimated as: 

d(l)y n ' \dQj \ a J T ' \EoJ J' ^ ' 



Pkink ■ Pph • Pmagn 
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An reasonable value for de^/dfjP' is about 0.1, that is for the fully ordered statg .the change in the dielectric function 
would be about 0.1. This is consistent with the ellipsometric data on NaV205LJ'EZI. Eq.(7) does not contain an evident 
small parameter which could allow to discriminate reliably intensities of the different contributions. Dominating 
phonon peaks arise due to their small damping F <C f2 ~ J rather than due to high efficiency of the phonon Raman 
scattering. For realistic parameters de^/dQ = 1, {zq/o)^ = 10~*, and F = 2 meV, pph at the phonon peak is of the 
same order of magnitude as pkink at i? ~ 2i?iat . This corresponds well to the experimental datallj, especially for the 
light polarized along the rungs, assuming that the continuum in NaV205 is formed by the kink-antikink excitations. 
The two-magnon mechanism is responsible for the Raman background in LiV2 05, where the kink-antikink pairs 
cannot be excited due to lone=range order of nonequivalent V"'+ and ions. In this compound the continuum 
is weak compared to phononsEj, however, their relative spectral densities are of the same order of magnitude as in 
NaV205 when the incident and scattered light are polarized along the ladders. 

At T > Tea the low-energy electronic excitations are long-ranged overdamped fluctuations of the order parameter 
coupled to the dynamical lattice distortion. These excitations have the same mechanism of Raman scattering as the 
kink-antikink pairs considered above, that is modulation of the polarizability due to intra-rung charge fluctuations 
leading to a broad intense continuum. The fluctuations of the order par-aHieter in the form of dynaniieal charge 
ordering persist up to T ~ 80K, as it was observed in Raman spectroscopyt2l, structural X-ray scatterin^J and spin 
resonanceE3 experiments. However, at present, it is hard to separate the contributions of the charge fluctuations and 
the two-magnon scattering to the Raman continuum. 

In conclusion, we investigated the low energy kink dynamics and Raman scattering by kink-antikink excitations 
in the charge ordered phase of a quarter- filled ladder with strong coupling to the lattice and considered NaV205 as 
an example. The spectral density of the scattered light and the overall scattering intensity by kink-antikink pairs 
very strongly depends on the charge order parameter. The mechanism corresponds well to three observed features 
of the continuum: (i) spectral range, (ii) depletion at Raman shifts less then 30 meV, and (iii) relative intensity to 
phonons. However, it is not possible at present to distinguish contributions of the kink-antikink and the conventional 
two-magnon processes. 
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FIG. 1. (a) Displacements of V ions (proportional to the length of arrow) in the zig-zag ground state and in a kink-like 
excitation. The radii of the grey circles correspond to the charges of the V ions, (b) The same for the in-line ordering, presented 
for comparison to (a). 




(a) (b) 

FIG. 2. The order parameter <;/>(?/, t) in kink-antikink pairs for different kink energies (upper panel) and absolute values of 
the changes in the dielectric function (Eq.(3)) caused by the pairs (lower panel) for time ti (solid line) and t2 > ti (dashed 
line), (a) Kink energy Ei, (b) kink energy E2 > E\. 
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